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Abstract 

We note that off the quark mass shell the operators (pi +Pf) tl j5 and ia fll ,(pi — 
PfYlhi both of which reduce to —a • E in the non-relativistic limit, are no 
£\j ■ longer identical. In this paper we explore the effects of this difference in the 

^o . 

contribution of these quark electric moments to hadronic electric moments. 
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I. INTRODUCTION 

When the loop diagrams which contribute to the electric dipole moments of quarks are 
evaluated in theories with P and T violating interactions, in general one finds a P and T 
violating quark-photon vertex function T^ q of the form 

f * = [D hq ( Pl + P/ )" 76 + iD 2>q a^{ Pi - p f ) vls ] Q q . (1) 

For convenience, we have extracted a factor Q q , the charge of the quark, in units of the 
positron charge. Heretofore it has been the standard practice to use the 75 form of the 
Gordon identity, valid on mass shell, 

uipf) [{Pi + PfT J5} u{pi) = u{p } ) [la^q^s] u( Pl ), (2) 

with 

<?=$-$, (3) 

to write the quark electric dipole moment interaction in its on mass shell form 

fv*.««~ = iD q Q q a^{ Pl - p f y l5 , (4) 

and to identify d q = D q Q q = [-Di, g + D 2 , q ] Q q as the electric dipole moment of the quark. 
Note that the dimension of d q is M _1 , and instead of introducing d q one could have made 
the natural units explicit and introduced a gyroelectric ratio h q , related to d q by d q = 
h q Q q e/(2m q ). However, unlike the practice in dealing with magnetic moments, it is standard 
practice to use d q rather than h q . For further convenience we also write di A = Di A Q q . 

In a previous paper [[j], the quark contribution to the electric and magnetic dipole mo- 
ments of the rho meson were calculated using vertices and propagators obtained from studies 
of the Dyson-Schwinger and Bethe-Salpeter equations. This provided some insight into the 
effect that QCD and confinement had on the electromagnetic properties of hadrons. The 
electric dipole moment of the quark used in that calculation was the on mass shell form of 
Eq. ([§. 



It should be noted that that the T- violating NN'y vertex of equation (||) is gauge invariant 

only when p\ = pi. However, this is enough for us to obtain a gauge invariant result for the 

P and T violating p — 7 coupling. As an off shell amplitude, it need not be gauge invariant 

- but it must satisfy the Ward Takahashi identities. We will examine the consequences of 

this requirement below. 

In this paper we extend that study to examine the effects of the vertex function ambiguity 
off shell, and calculate the electric dipole moment of the rho meson which is generated by the 
P and T odd vertex function T^ of Eq. ([!]). We continue to use the semi-phenomenological n- 
point functions of Quantum Chromodynamics obtained in the Dyson- Schwinger and Bethe- 
Salpeter framework. For ease of comparison with previous results we also give perturbative, 
bag model and non-relativistic results, as in our previous paper Q. The methodology 
employed here allows us to investigate to what extent, when quarks are confined within a 
hadron, the two terms give equivalent results. 

We emphasize that the rho electric dipole moment is studied as a model problem. The 
hadron electric dipole moment which has been extensively studied experimentally and theo- 
retically is the neutron electric dipole moment. The relevant calculations are reviewed in ||. 
We suggest that our results for the rho meson give an indication of the effects that one may 
expect for other hadron systems, including the neutron, when one has available techniques 
for exploring the off mass shell behaviour of the wavefunctionsQ. 

For a review of the topic of Dyson-Schwinger equations and their application to hadron 
phenomenology see ||. For examples of the use of the Dyson-Schwinger and Bethe-Salpeter 
equations in the calculation of hadronic properties see references within |],|3| . 

As in our previous paper, use will be made of forms for the quark propagator which 
ensure confinement by having no Lehmann representation, i.e. the quarks cannot go on 
mass-shell. Three different phenomenological forms for the rho amplitude will again be 



l K start was made on that investigation by Costella 



used. In that paper we used a dressed quark-photon vertex to ensure that both the Ward- 
Takahashi Identity and multiplicative renormalizability are preserved. In our present work 
we concentrate on the contribution of the P and T violating vertex function of Eq. ([[]), 
which we do not dress. 

II. THE INGREDIENTS OF THE CALCULATION 



A. The pp7 Vertex 

The form for the ppy vertex is given by JTIJ 



Vatf = d a »{p) [P^-guaSiq 2 ) + quq«Q{q 2 )) 

+{9^qa ~ 9^qu)M(q 2 ) + e vaw q p V(q 2 )] d/(p'), (5) 

with p the initial momentum of the rho, p' the final momentum and q the momentum of the 
photon. The d a p are the rho spin projection operators 






d a p{k) = g af3 5L_. (6) 



The form factors £ , Q, M. and T> are interpreted in the limit as q 2 — > through e£(0) = e, 
the charge of the rho, eJA(0)/(2m) = /x, the magnetic moment, Q(0) is related linearly to 
Q, the quadrupole moment, and eD(0)/(2m) = d, the electric dipole moment of the rho.[| 
The integral for the rho EDM is too difficult to perform analytically. This means that to 
isolate the relevant terms they will be projected out of the integral before it is integrated 
numerically. Using a projection operator given by 



2 We use m without subscripts for the mass of the p + . All other masses in this paper are identified 
with the appropriate particle by subscripts. 



P a » = -^dr f {pw a .df ia '{jf), (7) 



the following result holds 



P a ^V a ^^V(0). (8) 



B. Quark Contribution to the Electric Dipole Moment 

All calculations below are carried out in Euclidean space with co-ordinates xi,X2,xs,X4, 
and metric and gamma matrices given by g^ = 5^, •y fl = 7^ and {7^,7^} = 25^. Treating 
the u and (i-quarks as identical except for their charge, the impulse approximation to the 
quark contribution to the electric moment is given by 

f d 4 k 
4m/3 = ( _1 )y 7^J tr CFD[Tp(k;p - q)S(k- + )iT litq (k ++ ,k-+) 

xS(k ++ )T*(k + q/2 ] p)S{k+-)]. (9) 

where k a p = k + ^ + 4f, T^ refers to the rho meson amplitude and T^(k; p) = C^T p a (—k\ p)C 
with C the charge conjugation operator, 7274. t^ q is the P and T violating quark-photon 
vertex and S(p) is the dressed quark propagator for a quark of momentum p, all of which will 
be discussed below. The tr C FD operation is a trace over colour, flavour and Dirac indices. 
A tilde has been placed on I aii p to emphasize that this integral gives only the P and T 
violating rho-photon vertex. 

C. Quark Propagators 

The general form for the solution to the quark propagator Dyson- Schwinger equation B 



is 



S(p) = -% i>o v (p 2 ) + a s (p 2 ) 

= (ip'A(p 2 ) + B(p 2 )y 1 . (10) 



A model form for the propagator is given by ||[7| 

/•) _ „-bix\ /-) _ p -b 3 x\ j /i _ P -Ax~ 



b\X J \ b 3 x J \ V ^ x 
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x + fri 1 x 



m / 

7* -4- TT7.^ \ 



2(x + m 2 ) - 1 + e" 2 ^™ ) 
2(x + m 2 " 12 



ov(x) = — J ; - 2 x 2 mC™e ix , (11) 



where x = p 2 /2D, ay(x) = (2D)oy(p 2 ), 0"s = \/2Dos{p 2 ) and m = m/\/2D, D is a mass 
scale. (A = 10~ 4 is chosen to decouple the small and large spacelike p 2 behaviour in Eq. (|TID; 
i.e. to allow bo to govern the ultraviolet behaviour and b 2 the infrared.) The parameters 
C m , m, b ,...,b 3 are g], C m ^ = 0.0, C m=0 = 0.121, m = 0.00897, b = 0.131, 61 = 2.90, 
b 2 = 0.603 and 63 = 0.185, with the mass scale D = 0.160 GeV 2 chosen to give the correct 
value for /„.. Note that with these values of the parameters the mean light quark mass is 
5.0MeV. This form for the quark propagator is based upon studies of the Dyson-Schwinger 
equation for S(p) using a gluon propagator with an infrared singularity, 

g 2 D^(k) = (^ - y^j te*DS*{k), (12) 

and a dressed quark-gluon vertex MM. A sufficient condition for the lack of free quark 
production thresholds is the absence of timelike poles in the propagator. The model quark 
propagator given above is an entire function (except at timelike p 2 = —00) and so does not 
have a Lehmann representation. This means it can be interpreted as describing a confined 
particle and it ensures the lack of the unphysical singularities corresponding to free quarks 

in la^p- 

The electric dipole moment is also calculated using a form for the quark propagator 
developed by Mitchell and Tandy §J, to investigate p-u mixing. This propagator is given 
by 

a s {x) = C m e- 2x + ^ (l - e~ 2x ) , (13) 

av(x) = ^ ~£ 2 - 2X) - mC m e- 2 \ (14) 



To fix the parameters A = \JlD and Cm, a fit to (qq), f^, r n and the ir-n scattering lengths 
was done. With \{m u + ma) = 16 MeV a best fit was obtained for A = 0.889 GeV and 
Cm = 0.581 J9]. This form has a deficiency which can be seen in its failure to correctly 
model the behaviour of as away from x = 0, in the massless limit with dynamically broken 
chiral symmetry M. The large value for the mass of the quark used is related to this 



deficiency of the propagator ||T0"|j 



D. Rho Meson Amplitude 



The dominant Bethe-Salpeter amplitude for the rho meson is given by fTT| , [T2|| 






r'(M-7,-?r'i^, (15) 



where k is the relative momentum of the quark and anti-quark, p is the momentum of the 
rho meson and I and /i are flavour and Dirac indices respectively. This form ignores other 
allowable Dirac structure in the vector meson Bethe-Salpeter amplitude and so introduces 
errors of the order of 10% fT5| . Using the quark propagator defined in Eq. (|TT| ) and a 
Ball-Chiu quark-photon vertex (see later), Chappell uses the following approximate form 
brr,(* jf ) H 

r, = *-"'">' + y^V. (is) 

with cli = 0.38845GeV, a 2 = 0.01478, a = 2GeV. The values for the parameters were found 
by fitting to the experimental values of f p and g p7T7T ^i . 



Pichowsky and Lee, also using the quark propagator defined in Eq. ([TT|) and a quark- 
photon vertex of the Ball-Chiu type, use an identical form for T p , given by ]15 



r#, P )— »/4 + _|_, (17) 

where a v = 0.400GeV, b v = 0.008GeV and c v = 125.0. Note that the coefficient of the 
rational term in the vertex function, and the scale factor, differ dramatically in the two 
attempts to fit this form of vertex function, 



Mitchell and Tandy, with S(p) defined by Eqs. (|T3| and |14]), use a form for the amplitude 
T p given by 



T p (k,p) = e 



-k 2 /a 2 



(18) 



with a = 0.194 GeV. The momentum scale in the exponential has been reduced by a factor 
of 2 and the rational term in the vertex function has been omitted. 
The normalisation for the rho amplitude is fixed by |~5 

d A k dS{k^ 



/ x PaP/3 



N c tr D 



-r pa (k,p)s(k-)r p p(k,p) 



(2tt) 4 dp, 
+N c tr D J -^-S(k + )T pa (k,p)^^T p ,(k,p), 



(19) 



where k a = k + 2£. This condition, with the fact that the quark-photon vertex obeys the 
Ward Identity, ensures that £(q 2 = 0) = 1, i.e. that the rho has unit charge 0. 

E. The Quark-Photon Vertex 

The quark-photon vertex also satisfies its own Dyson- Schwinger equation, but solving 
this integral equation is difficult. Despite this, a realistic ansatz for the vertex function has 
been developed [pT5| — |TSf| - 

The quark-photon vertex ansatz thus obtained is given by 



^BC+CP/ 



ir +u *(p,q) = Q q r^(p,?) + r^(p, g ) 



(20) 



The Ball-Chiu vertex, r|^ c has the form [FJ 



rffeg) 



A{p 2 ) + A{q 2 ) (p + g)„ 

^ + ~2 



2 r 

,2\ of„2\ 



Q- 



-{A{p 2 ) - Aiq 2 ))^- 



- t(B(p 2 ) - B(q 2 )) 



(21) 



This vertex ansatz is completely described by the dressed quark propagator and satisfies 



both the Ward-Takahashi and Ward Identities, is free of kinematic singularities as q 2 



p 



transforms correctly under appropriate transformations and reduces to the perturbative 
result in the appropriate limit. 

To ensure multiplicative renormalizability Curtis and Pennington added a transverse 
piece to the Ball-Chiu vertex [[HJ. This term has the form 

r c P(Pi q) = i ~>^-^)-^,w-^ {A(f?) _ AW)) (22) 

with 

,2 



dip, q) 



p 2 + q 2 



(p 2 -g 2 ) 2 + (M 2 (p 2 ) + M 2 (g 2 )) Z j, (23) 

M(p 2 ) = ^1. 
yp ' A(p 2 ) 

We introduce the quark electric dipole moment through the vertex function T M of Eq. ([[]). 
We have not attempted to "dress" this vertex by including gluon loop corrections. 

To have an explicit example in mind we note that in models where the quark electric 
dipole moment is generated by Higgs exchange |2(]|n]], the structure of the quark electric 
dipole moment is that of Eq. ([!]). With the charged Higgs-quark interaction Lagrangian 
given by 



£ int = 2^Gfd 



V KM M D aH + R + MuV KM l3H + L D + H.C. , (24) 



the electric dipole moments of the up and down quarks are dominated by the bottom and 
the top quark loops respectively. With a charged Higgs mass of the order of the top mass, 
the electric dipole moment of the down quark is significantly larger than that of the up 
quark, and the P and T violating d — 7 vertex is given by Eq. (|l|), with 

d 14 = K (FxK/M 2 ,) + F 2 (m 2 t /M 2 H )) (25) 

d 24 = -KF l {m 2 t /M 2 H ) (26) 



K=^f\V dq \ 2 Im(a*P) (27) 



F 1 = i-^—(l + -^—hix) (28) 

1 — x \ 1 — x ) 

* = 7i^( fe - 3+? 7^rH' (29) 



Note that, in the on mass shell approximation, the F± terms cancel, leaving only the F2 
terms. 



F. Gauge Invariance, Ward Identities, Renormalization and Off-Shell behaviour 



We have already noted the the amplitude of ([!]) is not manifestly gauge invariant, since 



<^f> = dM - P }h 5 



(30) 



which vanishes only when p\ — pi, and thus vanishes on mass shell. However, because f £ is 
a Green's function, not an element of the S'-matrix, it need only satisfy the Ward-Takahashi 
identity |22|], 

<9X 5 (p) 



t m (p,p) 



dpu 



(31) 



or 



q^iPhPf) = ~ (S 5 (pj) - S 5 (P/)) • 



(32) 



The CP- violating term, E5Q9), in the quark self energy, required to satisfy the Ward Identity, 
will be generated by the same process that generates the T violating vertex, for example by 
Higgs loops. In general 



E 5 (p) =m 5 (p)75, 



(33) 



and E 5 (p) is divergent, coming from a one loop diagram. There are however no terms 
of this structure in the unrenormalised lagrangian, so we have no parameter to adjust in 
the renormalisation process. In QED this is not a problem, as one can demand that the 
equivalent term vanish on mass shell |23|| . In the present case quarks do not have a mass 
shell in the confined phase of QCD, and the best we can do is to renormalize S 5 such that 
it vanishes at some renormalization point, /i. That is, we have, after renormalisation, 



S 5 (p) = -di(p 2 -^ 2 )7 5 , 



10 



(34) 



guaranteeing the Ward-Takahashi identities^. For non-confined particles the choice \x = m, 
where m is the mass of the particle concerned, recovers the results of [^3| . 

To lowest order then, the corrections to the quark-photon vertex due to Higgs loops are 
given by the sum of diagrams given in Fig. ([[]). If p is set equal to m q then the "non- 
standard" quark-photon vertex term in Eq. ([I]) is cancelled, with the result that 

f£ = (d 2 + diV'-Vys, (35) 

just as if we had used the on-shell identity of Eq. (§). If, on the other hand, one renormalizes 
S5 at some other point p 7^ m q , then there is no such cancelation but T^ has a p? dependence. 
This p 2 dependence is a result of the absence of an adequate theory of perturbations about 
the confined phase, indeed of an absence of an adequate theory of the confined phase. This 
being the case, to make progress we set set p 2 = {p av ) 2 , with the result that the last two 
diagrams of Fig. ([I]) give negligible contributions and so 

f£ = di(p< + P/) M 75 + d 2 %a^q vlh , (36) 

just as Eq.([I|). 

It may be thought that the lack of gauge invariance of the CP violating vertex of ([!]) 
would lead to a lack of gauge invariance for the electric dipole ppj coupling. Explicit 
calculation shows this not to be the case -- the resulting electric dipole pp'-f coupling has 
only the gauge invariant form e vapbp q p . While a calculation of the neutron EDM using the 
T and P violating vertex of Eq. (|l|) in a three quark model of the nucleon is at present 
impractical, if one uses it in a quark-diquark model of the nucleon, the resulting NN^y, T 
and P violating, vertex is similarly gauge invariant. 



3 At least when we can regard d± as independent of pi and pf 
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III. THE DIPOLE MOMENT 

The dipole moment of the rho can now be calculated. The colour, flavour and Dirac 
traces are performed for the integral I a ^/3, see appendix A, and the dipole term projected 
out using the projection operator P a ^ . The integral is then performed numerically using 
Gaussian quadrature methods to obtain the results given in Table |. 

For comparison, we also use perturbation theory to calculate the electric dipole moment 
of a rho meson made of free, undressed quarks. To do this for the (pi +p/) M 7s termQ in the 
integral I af i/3, the following replacements are made, 

S(p) -► l/(ij + m) (37) 

r"(M-i (38) 

t' fl (k ++ , fc_ + ) -► (k ++ + k_ + )^ 5 D hq Q q . (39) 

Using the standard dimensional regularization scheme along with the Feynman parameteri- 
zation technique |23j , the integral can be shown to remain finite. However, as was shown in 



our earlier paper 0, the corresponding contribution to the electric charge is infinite. Thus 
calculating the contribution to the rho dipole moment in perturbation theory via the ratio 
£>(0)/£(0) yields 

d p;peIt = 0. (40) 

In QjJ the rho EDM was also calculated in the bag model. As it has been noted previously 
that ip(pi + Pf)(i,j5ip = ipio fiuQ 1 'l^ip f° r quarks on the mass-shell, it follows that the non- 
standard form for the quark EDM will yield the same bag model result as the standard 
term, since the quarks satisfy the equation of motion inside the bag. 



4 the equivalent calculation for the icr^ v q u ^ contribution was performed in our earlier paper [|J, 
with the result that the rho electric dipole moment is d p = {rriudu+mddd) / '(2m) in this perturbative 
limit. 
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In the non-relativistic limit, the result is simply 

d { p NR) = {(di, u + d 2>u ) + (d ltd + d 24 )} . (41) 

The results may be expressed in the form 

d p = (d ljU + d 14 ) A + (d 2 , u + d 2 , d ) B (42) 

For our various choices of input, A and B are given in Table |. 

IV. DISCUSSION 

Using a model form for the quark propagator obtained from the Dyson- Schwinger equa- 
tions, two different phenomenological rho-meson amplitudes fitted to f p and g p7Tn and one 
used to describe p-u mixing, and a quark-photon vertex which incorporated the P and T 
violating quark-photon vertex function of Eq. (Q), we have calculated the electric dipole 
moment of the rho. 

For each of the coefficients A and B, all three non-perturbative models agree qualitatively 
with the Bethe-Salpeter amplitude of [|1| yielding the largest value for d p and that of || 
the smallest. The results for the amplitudes of Chappell, and Pichowsky and Lee, are very 
close, which is indicative of the similarities of the two approaches. 

For the same non-perturbative model, a comparison of the coefficients A and B shows 
that the contribution of the (pi + P/) Ai 7s term is rather less than that of the "standard" 
^/ii/9^75 term. The (pi + P/) Ai 7s contribution was ~ 67% of the io pv (f^ result for the 
amplitude developed by Mitchell and Tandy, ~83% for that of Chappell and ~81% for the 
amplitude of Pichowsky and Lee, with the (pi+Pf)^ and ia pv q u ^ terms weighted equally. 

In the Bag Model, and in the non-relativistic limit, these two terms contribute equally. 
In the perturbative, free limit, on the other hand the ic pv q Vr )^ term gives a finite (but small) 
contribution, whereas the (pi + p/^75 term gives a vanishing contribution. In a sense the 
non-perturbative models interpolate between these limits. 

13 



To be explicit, consider the Higgs model introduced above. We take the Higgs mass 
to be lOOGeV, and the top mass to be 170GeV, and note that, in this model, d d ^> d u . 
Eqs. (|25|)-(|29|) show that the (p; +VfYlb term has a coefficient more than twice as large as 
that of the standard term, i.e. di^/d2,d = 2.34. The resulting rho electric dipole moment, 
expressed as a fraction of the non-relativistic value is given in Table [IT]. 

In this case the models of Chappell, and Pichowsky and Lee, give similar results which are 
about two thirds of the non-relativistic result. The Mitchell and Tandy result is significantly 
smaller. The good news is that the standard non-relativistic result is good to a factor of 2 
to 3, and that while one is concerned with approximate estimates it is a reasonable guide to 
the overall effect. 

Though this calculation was carried out for the rho meson it leads to some interesting 
speculation about the calculation of the quark contribution to the neutron electric dipole 
moment, which is of interest in the study of CP-violation. The use of the full off mass shell P 
and T violating interaction of Eq. ([I]) could significantly alter both the standard model and 
the many non-standard model results for the neutron EDM . While one can expect the order 
of magnitude of the existing calculations to be correct, at a precision level the off-mass 
shell effects considered here will become important. As some of the non-standard model 
calculations yield values that are close to the current experimental limits [@], it will become 
important to attempt precision calculations in the future. Of course, when such calculations 
are attempted, many additional sources of a neutron electric dipole moment will need to be 
considered. In the Higgs model, these include 

1. charged Higgs contributions to 

(a) the quark colour-electric dipole moment 

(b) the gluon colour-electric dipole moment 

(c) other P and T odd quark-quark interactions 

2. neutral Higgs contributions to the above 

14 



3. hadronic effects 

Considerations similar to those discussed here will enter the calculation of the quark colour- 
electric dipole moment effects, and possibly also the other additional effects. 

It is well known that the fact that the quarks in hadrons are both relativistic and bound, 
and thus are significantly off their mass shell, has significant effects on the quark contribution 
to the magnetic moments of hadrons |flj,|25|l. The study of the rho meson presented here 
shows that similar effects occur in the calculation of the electric dipole moments of hadrons. 
In particular, it is necessary to differentiate between the two forms of the P and T odd 
quark-photon interaction in these calculations, and not to treat them as being equivalent as 
they are on the quark mass shell. The next generation of calculations of the electric dipole 
moment of the neutron should take such effects into account. 

APPENDIX A: INTEGRAL FOR THE ELECTRIC DIPOLE MOMENT 

Taking the colour and flavour traces, the integral becomes 
2iN r d 4 k 

I** = ^rj j^r P ( k ) T p( k + <i/ 2 )tr D 

+7/3 ^-+75 #++7aVri + 7,3 /f?-+757a ^+-k^T 2 

+7/3 f-^l^ W (f H++laT 3 + 7/3 ^-+l^ W q P la ft+-T 4 

+7/375 #++7a ft+-k^T b 

+7/375<Vp? P /^++7a ft+-T 6 + l(3l5^^ P q P laT 7 , (Al) 

where k^ = (k ++ + A;_ + ) M , only P and T odd terms have been retained, and terms which 
will not survive the Dirac trace have been omitted. The Ti_ 7 are defined by 



Ti = -V l a v {k I _ + )a v {ki + )a s {K_), 
T 2 = -V^vikt^asikl^avikl.), 
T 3 = -V 2 a v (k 2 _ + )a v (k 2 ++ )a s (k 2 + _), 
T 4 = -V 2 a v (k 2 _ + )as(k 2 ++ )a v (kl_), 
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with, 



2 



T 5 = -V 1 a s (e + )a v (k 2 ++ )a v (kl_) 
T 6 = -V 2 a s (k 2 _ + )av(kl + )a v (kl_) 
T 7 = V 2 a s {kt + )a s {k 2 ++ )as{kl_), 



(A2) 



Vi = - (di,„ + di jd ) , 
K 2 = -i (d 2jU + rf 2 ,d) • 



(A3) 
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TABLE I. Expansion coefficients A and B for the Electric Dipole Moment d p 



Model 


A 


B 


Chappell 


0.620 


0.743 


Pichowsky and Lee 


0.627 


0.779 


Mitchell & Tandy 


0.418 


0.627 


Perturbative 


0.000 


0.010 


Bag Model 


0.828 


0.828 


Non-Relativistic 


1.000 


1.000 
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TABLE II. Electric Dipole Moment d p in a Higgs model, as a ratio to the Non-Relativistic result 



Model 


Ratio 


Chappell 


0.657 


Pichowsky and Lee 


0.672 


Mitchell & Tandy 


0.465 


Perturbative 


0.003 


Bag Model 


0.828 


Non-Relativistic 


1.00 
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+ 




+ 




FIG. 1. Lowest order Higgs loop corrections to the quark-photon vertex 
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